Jamming transition in traffic flow under the priority queuing protocol 
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Packet traffic in complex networks undergoes the jamming transition from free-flow to congested 
state as the number of packets in the system increases. Here we study such jamming transition when 
queues are operated by the priority queuing protocol and packets are guided by the dynamic routing 
protocol. We introduce a minimal model in which there are two types of packets distinguished by 
whether priority is assigned. Based on numerical simulations, we show that traffic is improved in 
the congested region under the priority queuing protocol, and it is worsened in the free-flow region. 
Also, we find that at the transition point, the waiting-time distribution follows a power law, and 
the power spectrum of traffic exhibits a crossover between two 1//" behaviors with exponent a ~ 1 
and 1 < Q < 2 in low and high frequency regime, respectively. This crossover is originated from a 
characteristic waiting time of packets in the queue. 
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Information packet transport via the Internet is an im- 
portant problem in complex systems from the perspec- 
tives of both theory and apphcation. Data packets cre- 
ated at certain nodes in the Internet travel to their desti- 
nations under transmission control protocols. During the 
journey, packets interact with other packets as they share 
a common hne or buffer in the network. Accordingly, sev- 
eral types of cohective behaviors can emerge in the form 
of self-similar traffic P, Q or chaos 0] . To enhance trans- 
port efficiency, one would hke to design an appropriate 
protocol to transport as many packets as possible with 
the fowcst cost. With this goal in mind, several packet 
transport models have been introduced. 

When a packet is sent from one node to another in a 
network, it is usually routed afong the shortest path; such 
a path is undoubtedly the best route when the number of 
packets in the network is relatively small. However, when 
many packets arc floating around in the network, traffic 
congestion can occur. This problem can be especially 
serious in scalc-frec (SF) networks since hubs are the 
bottlenecks of traffic flow. To resolve this congestion, 
many routing protocols have been proposed, including 
the hub avoidance protocol Q and the optimal routing 
protocol 0- These are static routing protocols, so that a 
path from one node to another is fixed regardless of the 
traffic level in the system at any given time. In contrast, 
there is a dynamic routing protocol that guides packets 
to alternative paths depending on the traffic on the path 
to each target [1] . 

When a packet travels through a node (router), it is 
temporarily stored in the buffer (queue) at the node. 
There can be many queuing protocols that control the 
order of packet transmission in the queue. The most 
common one is the 'first-in-first-out' (FIFO) protocol. 
Alternatively, the 'last-in-first-out' (LIFO) protocol can 
be used [l0l.[lll . The priority queue is a rather differ- 
ent protocol [l2||. Each packet is assigned a priority upon 
its birth. A packet with the highest priority is treated 
first in the queue, irrespective of its order of arrival. The 
diffusion process (isj under the priority queuing protocol 



in complex networks have been studied previously. 

In this paper, we study the jamming transition of 
packet transport on SF networks using the priority queu- 
ing protocol and the dynamic routing protocol by adapt- 
ing Dijkstra's algorithm [iJl- At each time step, every 
node creates a packet with probability q whose destina- 
tion is chosen randomly. These packets are assigned to be 
either with or without priority. The fraction of the pack- 
ets that are priority-assigned is /. In the queue, pack- 
ets with priority are delivered first, followed by packets 
without priority. Similar types of packets in the queue 
are treated following the FIFO protocol. Packets with 
priority may be regarded as paid packets when down- 
loaded from a certain web site. Note that when / = or 
/ = 1 , the priority queuing protocol reduces to the stan- 
dard FIFO protocol. We present a phase diagram for the 
free-fiow and congested phases in the parameter space 
((?,/), and we show that the priority queuing protocol is 
efficient when the system is congested. 

We simulate packet transport under the dynamic rules 
below on undirected binary scale-free networks generated 
by the static model [l^. For the network, the total num- 
ber of nodes is = 1, 000, the average degree of a node 
is (fc) ~ 4, and the degree exponent is 7 = 2.5. Each 
type of packet travels along the path that minimizes the 
quantity 



(1) 



where ix is the hopping distance along a path x be- 
tween nodes s, d, Qi{t) is the queue length at node i 
on the path x, and /i is a traffic-control parameter 
For packets with (without) priority, Qi{t) is regarded as 
the number of priority-assigned (both types of) packets 
that have accumulated in the queue at node i. Hence, 
the path minimizing L^^d can be the best choice to route 
the packet at time t, since the packet can circumvent 
congested nodes along its way. This path is determined 
using Dijkstra's algorithm in the simulation. Note that, 
when /i = 0, this routing protocol reduces to the shortest 
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FIG. 1: (Color online) The delivery fraction D for packets 
with priority versus the packet generation rate q for various 
h. The traffic-control parameter h; h — 0.0 (□), 0.1 (0)i 0.5 
(A) and 1.0 {(}). The priority fraction / is chosen as 0.9. 



path routing protocol. We choose other control factors 
as follows: The queue size is unlimited, and the process 
rate in each queue is one packet per time step. Thus, if 
more than one packet arrives at a node per unit time, 
then the queue length increases. The system is updated 
in parallel, meaning that all packets move simultaneously 
with the queue length information of the previous time 
step. Our simulations are performed for up to 10* time 
steps. 

To characterize the jamming transition, we define an 
order parameter for each type of packet as the delivery 
fraction Da , where a = p for packets with priority, a = n 
for packets without priority, and a = tot for all packets 
combined. is defined as 



Da = lim 



1 



t^ooNait-to) 



\a{t')dt', 



(2) 



where \a{t) is the number of packets of type a arriving 
at their destinations at time t and Na is the number of 
packets of type a generated in a unit time. Here, Na is 
Na = Nqf, Nq{l — /), and Nq for a = p, n, and tot, 
respectively. If the traffic of packets of type a is in the 
free- flow state, then Da is close to 1, since all packets 
of type a arrive successfully at their targets. However, 
if the traffic is in a congested state, then < Da < 1- 
When Da = 0, the traffic is in the completely congested 
state. 

In general, the jamming transition point for packets 
with priority does not coincide with that for packets with- 
out priority, denoted by qp and gn, respectively. Obvi- 
ously qp > qn- When f ~ and / = 1, all packets are 
of the same type and the queuing process is governed by 
the FIFO protocol. In this case, the delivery fraction and 
the jamming transition point are denoted by Dq and qo, 
respectively. 

First, we consider the effect of the traffic-control pa- 
rameter h. We observe that the traffic is dramatically im- 
proved when h > compared with the case when h ^ 0, 
as shown in Fig. [TJ However, the traffic seems to be un- 
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FIG. 2: (Color online) The delivery fraction D for packets 
with (a) and without priority (b), and all packets (c) versus 
the packet generation rate q for various fractions / of the 
priority assignment. The inset at (a) shows Dp collapses a 
unique scaling function. 



affected by changes in h when h > 0. This result is rea- 
sonable because the sum of accumulated packets along 
the path is far larger than the topological length in 
the congested state. Hence, we will confine our interest 
to the case when h = 1 from here on. 

We observe the behavior of the delivery fraction. We 
consider Dp as a function of q for various / in Fig. [2lja). 
The preliminary result of this was reported in [l6| . For 
/ = and / = 1, we obtain qo ~ 0.07, which can change 
depending on the system size N. When < / < 1, the 
jamming transition point for priority-assigned packets 
is larger than go- Since packets without priority do not 
hamper the traffic of packets with priority, one can obtain 
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FIG. 3: Phase diagram for packet traffic. Tlie dark grey 
region represents the free-fiow state when the FIFO protocol 
is used. The light grey region represents the additional free- 
flow state for the priority-assigned packets under the priority 
queuing protocol. 

the relation 

qpf = qo- (3) 

On the other hand, when q > q-p, a fraction of packets 
with priority cannot reach their targets during a given 
time interval. Hence, -Dp < 1 in such a case. 

Second, we examine the behavior of the jamming tran- 
sition for packets without priority. Since packets with 
priority delay the traffic of packets without priority un- 
der the priority queuing protocol, the jamming transition 
for packets without priority occurs at q^ smaller than go- 
As > qa for a given /, Dp = 1 and < < 1 
for gn < g < Qp. Fig. [2Ub) shows the behavior of £)„ 
as a function of q for various /. For q > qp, the traffic 
for packets without priority is completely congested, i.e., 

Next, we combine the above two cases and consider the 
jamming transition for all packets irrespective of priority 
assignment. Dtot is the order parameter for all packets. 
The behavior of Dtot is shown in Fig. [D^c). Dtot satisfies 
the relationship 

Act = fDp + (1 - /)A, (4) 

where Dp and can change depending on q. We sum- 
marize the delivery fraction for all packets as 

{1 for g < qn, 

/ + (1 - /)£>„ for qn<q< qp, (5) 
fDp for q> qp. 

It is interesting to note that the delivery fraction Z?tot 
under the priority queuing protocol can exceed the Dq ob- 
tained from the simple FIFO protocol. This phenomenon 
can occur in the region of g > qp, as shown in Fig. [^Jc) 
(for example, q > 0.12 for / = 0.5, and q > 0.096 for 
/ = 0.75). The unexpected improvement in overall trans- 
port efficiency is due to the priority queuing protocol, 
since it enables us to control the density of packets with 




FIG. 4: (Color online) The waiting-time distribution in the 
free-flow region for packets with (a) and without (b) priority 
for / = 0.9 and / — 0.25 respectively. The power-law behav- 
iors are obtained at qp — 0.053 and g,i = 0.015, which are 
regarded as the jamming transition points. The solid lines 
are guidance with slopes —3.4 (a) and —2.3 (b), respectively. 



priority; these packets remain deliverable packets in the 
congested state. For the other cases, Dtot < Dq, im- 
plying that the overall traffic under the priority queuing 
protocol is worse in the free-flow state, despite the im- 
provements observed when the system is in the congested 
state. 

Fig. [3] is a phase diagram of the traffic of packets in the 
space of (q, /). Under the FIFO protocol only, the phase 
space is divided into two parts by the line q = qo, with 
the free-flow state appearing in the region q < qa and the 
congested state appearing in the region q > q^. However, 
when the priority queuing protocol is used, the free-flow 
region can be extended into the region q < {qa/ f). 

We measure the waiting-time distributions P{t^) for 
packets with and without priority [l3| . The waiting 
time ty^ is defined as the time spent in the queues on the 
way to the target, excluding the transit time. As shown 
in Fig. m the waiting-time distribution follows a power 
law P{tvj) ~ near the transition points qp and q^, 
where S « 3.4(1) for packets with priority and S w 2.3(1) 
for packets without priority. In the free-flow regions q < 
qp and q < q^, the waiting-time distribution behaves as 
P(tw) e"*"/"^, where r is the mean waiting time and 
depends on the density of packets in the network. The 
maximum waiting time of packets in the free-flow region 
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FIG. 5: (Color online) The power spectrum of traffic in the 
system for packets with (a) and without (b) priority, re- 
spectively. We use / = 0.9 (a) and / = 0.25 (b). For 
= 0.053 (a) and = 0.015 (b), the power spectrum ex- 
hibits a crossover between two power-law behaviors. Solid 
lines are guidelines with slopes —1 in small uj regime and 

— 1.5 in large uj regime for (a) and —1 in small uo regime and 

— 1.8 in large ui regime for (b). 



is estimated to be t,„ « 70 for a system of size N ~ 10"^. 
This means that packets can reach their targets within 
tm at the most. 

The power spectrum S{ijj) of the traffic is defined as 



(6) 



where g{uj) = Yiit=o ^(0^ '"t'^* and F{t) is the number 



of packets with (without) priority in the system at time 
t. The power spectrum of the traffic in the system is 
measured in Fig. [H The behavior of the power spectrum 
S{ijj) depends on the packet generation rate q as well 
as the packet type. When the packet generation rate 
q is near the jamming transition point, the power spec- 
trum exhibits a crossover between two power-law behav- 
iors S{u)) ~ w"'' with w 1 and rj k, 1.5(1) for packets 
with priority and ry « 1 and -q w 1.8(1) for packets with- 
out priority. Such behaviors are different from what is 
observed in model systems and empirical data [l3, [13 ■ 
This difference is probably due to the short-tailed be- 
havior of P(tw) compared with that for [Tol . \v^ . This 
difference of P{tvj) is caused by the difference in rout- 
ing and queuing protocol. The crossover behavior occurs 
roughly at uj^ « T/2-nUn w 10^, where T = 5 x lO'^ is 
the total simulation time step. This value corresponds 
to the maximum waiting time in the system, tm, roughly 
estimated in Fig. 0] to be 60 ^ 70 for a system of size 
N = lO'^. The l//-type power spectral density suggests 
that there exists a long time correlation in the transport 
of both types of packets. 

We also perform the same simulations on the Erdos 
and Renyi (ER) network to ascertain whether our result 
is affected by network structure. With the exception of 
the increment of the transition point, the generic features 
of the simulation results remain unchanged using the ER 
network, indicating that our results are independent of 
network structure. 

In summary, we studied the packet transport problem 
on SF networks under the priority queuing protocol and 
the dynamic routing protocol. We showed that total traf- 
fic can be improved in the congested state by introducing 
the priority queuing protocol, although the overall traffic 
is worse in the free-flow state and the jamming transition 
point is reduced. The jamming transition points for pack- 
ets with and without priority are different. Near each 
jamming transition point, the waiting-time distribution 
follows a power law, and the power spectrum exhibits a 
crossover between two power-law behaviors. We obtain 
l//-type power spectra in the small lo regime for both 
types of packets. 
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